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RESEARCH MEMORANDUM 


COMPARISON BETWEEN THEORETICAL AND EXPERIMENTAL RATES 
OF ROLL OF TWO MODELS WITH FLEXIBLE RECTANGULAR 
WINGS AT SUPERSONIC SPEEDS 


By John M. Hedgepeth and Robert J. Kell 
SUMMARY 


A comparison is presented between the experimentally measured and 
theoretically calculated (by the method of NACA TN 3067) rates of roll 
of two rocket-propelled models with flexible rectangular wings. The 
comparisons show that although there are large aeroelastic losses in 
rolling rate, the theory predicts the actual rate of roll accurately. 
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INTRODUCTION 


In reference 1, a method is presented for calculating the aero- 
elastic effects at supersonic speeds on the rolling behavior of air- 
craft with flexible rectangular wings. The method employs structural 
influence coefficients to determine the deformations of the wings and 
linearized supersonic lifting-surface theory to find the airloads. 


The purpose of the present paper is to assess the accuracy of the 
method of reference 1 by comparing its predictions with some experimental 
data presented in reference 2 for two, rocket-powered test models. Com- 
parisons between theory and experiment are given in the form of plots of 
rolling rate against Mach number. 


SYMBOLS 
D local flexural stiffness, Et//12 (Gg 2) 
E Young's modulus of elasticity 
G shear modulus of elasticity 
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GJ elementary torsional stiffness, G dx 
-c/2 Í 
IGm(Y 52) 
——= structural rate-of-twist influence function which results 
oy from a unit concentrated torque 

I moment of inertia of beam that represents flange effect 
of aileron 

K(x) stiffness of root springs 

L(y) aerodynamic load per unit span, positive upward 

M free-stream Mach number 

M(y) aerodynamic moment, per unit span, about midchorã, 
positive in positive twist direction 

Ph static pressure at altitude 

Po standard static pressure at sea level 

Q(y) aerodynamic moment, per unit span, about elastic axis, 
positive in positive twist direction 

V free-stream velocity 

W(y) deflection of midchord line of wing, positive upwards 

a ratio between fuselage radius and exposed wing semispan 

b total wing span, 2(al + 2) 

c wing chord 

Ca aileron chord 

e distance measured forward from midchord to elastic axis, 
expressed as fraction of chord 

l exposed wing semispan 

P rolling velocity 

pb r 


— tangent of wing-tip helix angle, positive for counter- 
2V clockwise roll when viewed from behind 
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8 ratio of root flexibility to wing flexibility 

t thickness of wing cross section 

w(x,y) local deflection of wing, positive upwards 

X,J coordinate system 

x distance from midchorã to location of beam that represents 
flange effect of aileron 

9 (y) angle of twist of wing, leading edge up 

E twist at mid exposed span which results from a unit 
torque at the tip 

5 aileron deflection, positive down 

A = ‘nondimensional parameter -involved in the rate-of-twist 

influence function (see eq. (A10)) 
H ' Poisson's ratio | 
p | rolling effectiveness, (pb/2V)p 
(pb/2V)p 

Subscripts: 

F flexible wing 

R rigid wing 

rev É aileron Teyersa 


DESCRIPTION OF THE MODELS 


The two models investigated (the last two in table I of ref. 2) were 
essentially the same in size and shape, the only important difference 
being that one had aluminum wings and the other had steel. Both models 
had three rectangular wings equally spaced around the rear portion of a 
long cylindrical body. (See fig. 1.) The wings were uniform in the , 
spanwise direction and had NACA 65A00% airfoil sections. The full-span 
trailing-edge ailerons were formed by bending the wing along the 
80-percent-chord line. 
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The pertinent dimensions of the two models are given in table I. 
Included in table I are the experimentally determined values.of 6,, 


the twist at the midspan due to a unit torque at the tip. These values, 
which were obtained from reference 2, are used as an aid in determining 
the structural characteristics of the wings. The ordinates for the NACA 
65A003 airfoil, obtained from reference 4, are also included in this 
table. 


THEORETICAL RESULTS 


The experimental data in reference 2 are given in the form of plots 
of rolling rate pb/2V against Mach number M. Also given are the vari- 
ations with Mach number of altitude (specified in the form of static 
pressure) during the flights. It is desired to calculate theoretically 
the variation of pb/2V with Mach number for the prescribed altitude 
variation for each of the two models. In order to do this, the compu- 
tational procedure outlined in reference 1 is followed. 


In reference l, several alternative calculation schemes were 
described. The particular one used herein is the same as that used in 
the analysis of the example configuration in reference 1. This approach 
is exemplified by the matrix equations (29) or (41) of that report. 
These equations are written in terms of the rate of twist do/dy rather 
than the twist 6 itself and make use of an interpolation procedure to 
reduce the number of degrees of freedom involved. For the present prob- 
lem the structural ingredients of these matrix equations - the rate-of- 
twist influence functions - are derived in the appendix of this report 
by the application of an approximate plate theory and are tabulated for 
both models in table II; the aerodynamic ingredients - the various 
indicial loads and the loads due to roll and aileron deflection - are 
obtained from reference 1. (The assumptioris made that the loads on 
each wing of the three-winged aircraft considered in this paper are the 
same as those resulting from the two-winged configuration considered in 
ref. 1.) It should be noted that the analysis in the appendix indicates 
the existence of an “elastic axis," a line along which loads can be 
placed without producing any appreciable twist. Accordingly, equa- 
tions (29) and (31) of reference 1 are modified as suggested therein to 
take advantage of this elastic axis. 


By using these modified equations, then, and by following the sug- 
gested computational procedure, the theoretical results presented herein 
were obtained. The results for aileron reversal appear in terms 
of (Ph [Po) , the ratio between the static pressure at which the ailerons 

rev 


reverse and the standard sea-level static pressure. In figure 2 is shown 
the variation with Mach number of this ratio for both the aluminum and 
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steel wings. A calculation of the rolling effectiveness q = BF / (Ee) 
F R 


for other values of static pressure at several different Mach numbers 
showed that the variation of q with Py/Po was almost exactly linear. 


Therefore, the values of ọ for any pressure ratio at a given Mach num- 
ber can be deduced from the value of (Ph /Po) for that Mach number. 
rev 


Thus, for either model 





(Ph / Po) 


(Pn/ Po) rev 


The rolling rate pb/ 2V for the flexible wing can be obtained by 
multiplying q by (pb/2V)p. Consequently, 
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The variation of (g 3) with Mach number has been found by the 
r/R 
method of reference 1; these values were used in conjunction with the 
information in figure 2 and the plots of the actual flight values 
of Pa P b against Mach number from reference 2 to obtain the theoretical 


p=1- (1) 


values of Ble shown in figure 3 for the two models. Also shown in 
figure 4 is the theoretical a © for the rigid wing. It should be 


remarked that measured values of Ph Po were lacking for Mach numbers 


higher than 1.4 for the steel wing and 1.8 for the aluminum wing; accord- 
ingly, the curves for the flexible wings have been stopped at these values. 


COMPARISONS AND DISCUSSION 


The variation of experimental = 5 with Mach number, as obtained 
from reference 2, is shown in figure 3 for comparison. The estimated 
experimental rigid rate of roll, obtained by extrapolation from the 
flexible data by assuming a linear variation of pb/2V with the param- 
eter Pp: is also shown in this figure. 
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From figure 5 it can be seen that although there is a large loss 
of rolling effectiveness due to aeroelasticity, the theory does a good 
job of predicting the actual flexible rate of roll. 


CONCLUSION 


A comparison between theoretically and experimentally determined 
rates of roll for two: rocket-propelled models with flexible rectangular 
wings shows that the method of NACA TN 3067 is capable of yielding 
accurate predictions for the aeroelastic effects on the roll of super- 
sonic aircraft with rectangular wings. 


Langley Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., June 11, 1954. 
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APPENDIX 


DERIVATION OF TORSIONAL INFLUENCE COEFFICTENTS 


In order to solve the aeroelastic rolling problem it is necessary 
to compute the torsional influence coefficients for the wing. These 
quantities are derived in this appendix by means of the same approach 
as that of reference 4; that is, the deflections are assumed to be lin- 
ear in the chordwise direction and the principle of minimum potential 
energy is employed. 


The structure under consideration is shown in figure (a). It con- 
sists of a solid plate that is uniform in the spanwise direction with a 
bent-up aileron and a more-or-less-rigid attachment to the model body. 
Two factors prevent the analysis of this structure directly by the method 
of reference 4: (1) The bend along the aileron hinge line produces a 
flange effect so that flat-plate theory cannot be used; (2) the root of 
the wing cannot be considered to be perfectly clamped. The analysis is 
therefore performed for the equivalent structure shown in figure }(b). 

In this figure, the aileron has been unbent and the flange effect has 
been represented by a beam. The moment of inertia of the beam is assumed 
to be equal to the difference between the moments of inertia of the bent 
and unbent aileron and the beam is located at the centroid of the 
difference-in-moment-of-inertia distribution. Thus, 


I = tano | i (x) (x - 5 + ca)? dx (A1) 
re t(x) (x - =+ ca)? xdx 
Cae, 
x=-£ (A2) 


Also, in order to represent the effect of incomplete root clamping, the 
plate-beam combination is assumed to be mounted on springs which prevent 
displacement of the root but permit non-zero slopes in the spanwise 
direction. In the analysis to follow, the shape of the stiffness 
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distribution of these springs, which, for the present, is indeterminate, 
is assumed to be such that the resulting equations exhibit their sin- 

plest possible form; the absolute magnitude of the spring stiffness for 
each of the two models is then selected so that the theoretical twist at 


the midspan due to a unit torque at the tip matches the experimental value 
in table I. 


The potential energy of this equivalent structure subjected to the 
distributed lateral load p(x,y) is 


c/ 2 2 
ailea (9) sS 


c/2 2 1 nef2 
L x OW x = X,y)w(x 
a K( ) E ( 9) dx JA f Je p( Vv) ( sy) dxdy 


-C 
(a3) 
EL) 
where w(x,y) is the deflection, D(x) = ——““—— is the plate stiff- 
ê, 
12\l -u 
ness, end K(x) is the stiffness distribution of the root springs. 
In accordance with the assumption of linear chordwise deformations, 
Let 


w(x,y) = W(y) - xo(y) (Al) 


Then, upon performing the integration with respect to x, 
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x = vi R + ET) (Ww')e - 2 (ap + EIx)w"o" + 
O 

(a, + ERÊ) (o) + 20 - u) a, (0" 2 ay + 


= [a (11º 2 = 2k,W'e' + (0 Er 


l 
J, [E W(y) + M(y)o(y)] ay (a5) 
where 
an 7 x2-l p(x) dx 
-c/2 
c/2 
= x271 
kn A K(x) dx 
and 


are, respectively, the section lift and moment about the midchord. 


Minimization of the potential energy yields the following differ- 
ential equations, 
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(a, + EI) WEY - (a, + Erx)elV = L 
(A6) 
= (ag + BIx)wlY + (a: + Bre) g iV _ e(1 - u)ajo" = M 


and boundary conditions, 


ee EI)W"(O) - (ap + ETx)e"(0) - kyW'(0) + kp0!(0) = O 
(22 + EIX)W"(O) - (as + ERÊ) 0"(0) - koW' (0) + kx0'(0) = O 
(a1 + EI)W"(2) - (ag + EIx)e"(2) = 0 
(ag + BTz)W"(2) - lz + EI) o"(1) =0 


(21 + EL)W"'(2) - (ap + EIX)6"'(2) = 


| 
O 


(22 + EIX W"'(2) - (a; + ERÊ 9"'(2) + 2(1 - H)jmo'(t) = 0. 


For purposes of aeroelastic calculations » Only the twist 0 is 
important; therefore, it is desirable to eliminate W from equations (A6) 
and the accompanying boundary conditions. In order to do this easily, 
it is convenient to take advantage of the freedom of choice of the shape 
of the root spring stiffness distribution and assume that 


ky B Ko E k 
aq + EL as + EIX az + EIX 


Since, as has been mentioned, the overall spring stiffness is to be 
selected by duplicating experimental values of twist, this limitation 
on the distribution should have negligible influence on the desired 
aeroelastic results. 
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After eliminating W, the resulting equations are 





2 
as + 22 - Set PPL L om yao” =m a LEE, (AT) 


a + EL a, + EL 


(0) = 0 


9"(0) = —Ł— e'(0) 
a, + EI 


o"(Z) = O 


2 
az + EIXº - terah o - 2(1 - p)aje'(2) = 0 


ay + E 


Integrating equation (A7) once and rewriting yields 
2 l 
err Mora A | a(n) an (A8) 
1° GJ 


where the remaining boundary conditions are 


o(0) = O 


e"(0) 


i 


= o '(0) (a9) 


o"(2) 


0 
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and where 
12 2(1 - p)a,2® (a10) 
az + E - a ELE)" 
21 + EL 
GJ = 2(1 - Bay (A11) 
Q(y) = M(y) - ecL(y) (412) 
i as + EIX 


-Teer + ET) E 


a + EL 
e da 


k47 aN) 


The parameter GJ can be written 


12 B(x) 
GJ = oJ m 


This quantity then is merely the elementary torsional sitffness for 
a thin cross section. The quantity Q is the total moment about an 
"elastic axis" located e chords ahead of the midchord. The term 
"elastic axis" is justified in this case because, as can be seen from 
equation (A8), loads applied at x = -ec anywhere along the span would 
produce no twist. The parameter s expresses the ratio between the 
overall stiffnesses of the wing and of the root springs. 


In view of the existence of an elastic axis for this configuration, 
only the influence functions due to torque need be obtained. More spe- 


cifically, it is desired to calculate the rate of twist = at sta- 


tion y caused by the application of a unit torque at station mn. This 
Ò 

quantity, designated (yn) , can be obtained by solving equation (A8) 
NA 


with the appropriate substi iia ee a i |. and is given by 


CONFIDENTIAL 
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OG, 1 cosh a(i - g) 7 
oy yan) ~ GJ cosh À + SA sinh À (- 1 + cosh 02s 

SA sinh Ad) (y > nm) 
ÔGm 1 1 
ny eee a a E h - X) - 
Oy yon) GJ cosh À + BA sinh À pee ae ter a ( l 


L 


sheimh A = simh A @ - 2) om nf (y< n) 


(A15) 


All the quantities in these equations except s can be found 
directly from the geometry and material properties of the wing. The 
parameter s can be evaluated by equating the theoretical and experi- 
mental values of the twist at midspan due to a torque at the tip. Such 
a process yields the following formula: 


sinh A ( - 3) sinh A Š + 





EE ig | 
a= 2030.) A 
neta - - r 


It should be noted that by determining s in this manner, the 
torsional characteristics of the wing are duplicated cldsely. If the 
value of s were taken to be zero, the root would be completely clamped 
and the wing would be too stiff. If, on the other hand, a value of s 
equal to infinity were used, the root would be completely free to warp 
and the wing would be too flexible. For the aluminum wing, for example, 
taking s to be zero and infinity would yield, respectively, 6, equal 
0.956 x 107+ ana 1.262 x 107+ radians in.-lb, whereas the actual experi- 


mental value is 1.175 radians in.-lb. From these values, it can be seen 
that appreciable error could result from an improper selection of s. 


In order to evaluate the rate-of-twist influence functions, the 
quantities I and x (eqs. (Al) and (A2)) were computed for the nominal 
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aileron deflection of 5°; in this calculation, the thickness +(x) was 
assumed to vary linearly from its value at the hinge line (see table I) 
to zero at the trailing edge. The values of a, as, and az (defined 


following eq. (A5)) were computed from the airfoil profile data in 

table I; in this computation, the necessary integrations were performed 
numerically by using Simpson's rule. The quantities A, GJ, and e 

were then evaluated. (See eqs. (A10), (ALl), and (A13).) Finally, equa- 
tion (A16) was used in conjunction with the experimentally determined 
values of 6, to find s. A tabulation of the pertinent parameters for 


each of the two models follows: 





Steel Aluminum 
E a ee, es. ee ee al. we Ag hs as RT in ee a OO 0.0485 
Gis bei so ar Bla we ee ee Cw ee we E500 42,100 
Da tae icine SM, O eds. ah Bir GRP Ke OSs Be Sie S wee ee | ITO 8.3138 
Boks, ‘ey is Bose Aes oe Re Bh OR OS ae ee a OOL 0.289 


The values of OG.) dy for the aluminum and steel wings have been 
computed for O < Z< 1 in steps of 0.2 and O <i<1 in steps of 0.1. 


The results are given in table II in matrix form. These matrixes are 
used directly for the computation of the aeroelastic effects on roll by 
the method of reference Ll. 
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TABLE T.- DESCRIPTION OF MODELS 


(a) Basic information 
model 


10.61 
7.07 
0.236 


0.2 


10.6 x 10º 

0.555 

Ər, radian/in.-lb .. 0.392 x 10-4 81.175 x 10-4 
DD 4.9 





“Two values of 6, for the aluminum model were given in 
reference 2. The value herein is the correct one. 


(b) NACA 654003 airfoil ordinates 


.500 | .0001} 
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O 0.115132 
O 0.021085 
O 0.003865 
O 0.000731 
o) 


0.000259 


O 0.398616 
O 0.144975 
O 0.027491 
O 0.005219 
O 0,00102} 


O 0.000375 


TABLE IT.- RATE-OF-TWIST INFLUENCE-FUNCTION MATRIXES 


0.351742 
0.398052 
0.072897 
0.01336) 
0.002526 
0.000895 


0.572192 
0.436856 
0.082838 
0.015727 
0.003085 
0.001129 


(a) Steel Wing, 


0.396851 
0.682758 
0.195062 
0.035760 
0.006760 
0.002396 


0.416156 
0.804593 
0.481008 
0.088182 
0.016671 
0.005908 


GJ = 117,500 1b-in.º 


O.4oh417 
0.856754 
0.767011 
0.210716 
0.039837 
0.014117 


0.429471 
0.888629 
0.941967 
0.783847 


0.220480 


0.078132 


(b) Alumimm Wing, GJ = 42,100 lb-in.* 


0.647775 
0.725521 
0.206553 
0.059215 


0.007693 
0.002816 


0.680688 
0.848060 
0.489192 
0.092874 
0.018220 


0.006670 


0.695021 
0.902382 
O. 771643 
0.215979 
0.042370 
0.015510 


0.703991 
0.936378 
0.948408 
0.781750 
0.225155 
0.082422 


0.430127 
0.892773 
0, 964701 
0.908006 
0.516103 
0.182891 


O. 705195 
0.940933 
0.972092 
0.906593 
0.517578 
0.189469 


0.430426 
0.694659 
0.97 504% 
0.964492 
0.918371 
0.427698 


0.705750 
0.94307 

0.983085 
0.964537 
0.823168 
0.435303 


0.430595 
0.895727 
0.980904, 
0.996499 
0.999338 
0.999765 


0.706073 
0.944.267 
0.989432 
0.99799 
0.999606 
0.999856 
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Figure 1.- Test model. 
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Figure 2.- Variation of pressure ratio at reversal with Mach number for 


ahe ue TWO noA Sigga 









20 NACA RM L54F23 





Theory 


Experiment 





Mach number, M 


Figure 5.- Variation with Mach number of the theoretical and experimental 
rate of roll. 
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(a) Wing with bent-up aileron and flexible root attachment. 
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(b) Equivalent plate-beam combination mounted on root springs. 
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Figure 4.- Wing with bent-up aileron and equivalent plate-beam configuration. 
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